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MEDIUM 

Yangsen Yao 1 , Norbert Schulz 1 , Q. Daniel Wang 2 , and Michael Nowak 1 

The Astrophysical Journal, 653, L000-L000, 2006 December 20 

ABSTRACT 

The iron depletion level and the gas-phase iron abundance in the hot (~ 10 6 K) interstellar medium (ISM) 
are critical to our understanding of its energy balance as well as the thermal sputtering, cooling, and heating 
processes of dust grains. Here we report on the first detection of the Fe xvn absorption line at 15.02 A from 
the hot ISM in the spectrum of the low mass X-ray binary 4U 1820-303 observed with the Chandra X-Ray 
Observatory. By jointly analyzing this absorption line with those from O vn, O vm, and Ne ix ions in the same 
spectrum, we obtain an abundance ratio as Fe/Ne=0. 8(0.4, 2.1) in units of the Anders & Grevesse solar value 
(90% confidence intervals). We find that the result is robust with respect to different assumed gas temperature 
distributions. The obtained Fe/Ne abundance ratio, albeit with large uncertainties, is consistent with the solar 
value, indicating that there is very little or no iron depleted into dust grains, i.e., most of or all of the dust grains 
have been destroyed in the hot ISM. 

Subject headings: ISM: abundances — ISM: dust — X-rays: ISM — X-rays: individual (4U 1820-303) 



1. INTRODUCTION 

As one of the most abundant refractory metals, iron is an 
important constituent of the interstellar dust grains and is be- 
lieved to ha ve the greatest frac tion of its atoms depleted into 
dust grains (lSofiaetal.l [r994). Surveys in ultraviolet (UV) 
wavelength band indicate that > 70% of the iron in the cool 
and warm medium of the Galactic disk/halo could be locked 
up in solid grains, and only the remains can be probed through 
absorption lines of gaseous Fe n (see Savage & Sembach 1996 
and references therein). Recent studies on X-ray absorption 
edges also provide evidence of iron depletion into dust grains 
in the interstellar medium (ISM; Juett et al. 2006). 

Iron can be liberated from the dust to the gas phase 
through thermal sputtering. Theoretical calculations show 
that the sputtering caused by high velocity shocks (v s ~ 50- 
200 kms" 1 ) could destruct several xl0% of grains in the 
ISM, and that with the presence of a magnetic field, grain- 
grain collisions can also be a very efficient process for dust 
grain destruction, specifically in the case of low-velocity 
shocks (vj<100 kms" 1 ; Draine & Salpeter 1979; Jones et 
al. 1994). The dust grains can also be fragmented/destroyed 
through thermal evaporating in the hot ISM. The recycling 
of atoms back to the gas phase has been evidenced by the 
difference of the elemental abundance ratios in the Galac- 
tic halo from those in the Galactic disk, and by the ioniza- 
tion dispar ity of the shocked material a t various post-shock 
distan ce dJenkins & Wallerstein 1 119951: ISavage & Sembachl 
1996). A comparison between the observed abundance ra- 
tio of (Mg+Fe)/Si in the dust of the halo clouds and that of the 
theoretical expectation, together with the observed correlation 
between the gas-to-dust ratio and the dust mass carried by Fe, 
indicates that the dust grain cores likely contain iron oxides 
and/or metallic iron. Some of these resilient cores can survive 
from t he shock destruction during the dust processing in the 
ISM (ISembach & Savage lfl996l[Frisch & Slavin II2003I) . 
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In the hot ISM, the knowledge of the exact amount of iron 
contained in dust grains, which is reflected in the gas-phase 
abundance of iron (e.g., Savage & Sembach 1996), is of great 
importance for understanding many astrophysical processes. 
For an emission spectrum of a solar abundance plasma at tem- 
perature ~ 10 6 K, the contributions of iron ions encompass 
about 50% of the total emitting energy and photons in the en- 
ergy range from 10 eV to 2 keV. On the other hand, dust grains 
are an effective coolant of the hot gas. The thermal energy of 
the hot gas can be transferred to dust grains via the collisions 
between electrons and dust grains, resulting in bulk heating of 
the dust and infrared (IR) emission. For a solar abundance of 
the dust grains, the grain IR radiation is at least comparable to, 
and could be more than 10 times more efficient than, the gas 
X- ray emission in cooling the hot gas at temperature ~ 10 6 
K dDwek & Arendtlll992l) . Consequently, the existence of 
the iron-bearing dust in the hot ISM could alter the chemical 
composition of the hot gas, therefore, largely change the flux 
and the spectral shape of the hot gas radiation, affect the en- 
ergy balance, change the cooling rate thus the lifetime of the 
hot gas, and eventually, adjust the star formation rate, in the 
whole ISM and affect the galactic evolution in general. 

The iron depletion level in the hot gas can be measured 
by modeling the iron emission lines or the absorption lines 
that the gas-phase iron ions imprint on the background point 
source spectrum, then comparing the iron abundance with the 
solar values. Directly measuring the IR emission of the dust 
in such environments is difficult due to the confusion with 
foreground cold dust in the Galactic disk along the line of 
sight. A recent study of the very soft (0.25 keV) X-ray diffuse 
background emission in the vicinity of the Su n suggests that 
the ga s phase iron is ~ 30% of the solar value (ISanders et al. I 
2001). However, given that the Sun may reside in a privileged 
location (i.e., the Local Hot Bubble), this depletion level may 
not be typical for the general hot ISM. Furthermore, the in- 
ferred emission measure is also sensitive to different adopted 
plasma emission models, and the less well known interac- 
tion between the solar wind ions with the local neutral ISM 
may further complicate the interpretations (e.g., Sanders et 
al. 2001; Pepino et al. 2004; Hurwitz et al. 2005). Absorption 
line studies, on the other hand, measuring the total column 
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TABLE 1 

Strong Transition Lines of Fe xvii and Fe xviii 
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FIG. 1. — Ionization fractions of oxygen, neon, and iron ions as a func- 
tion of temperature for a g as in the collisional ionization equilibrium state 
iArnaud & Ravmond 1992). The fractions of neon and iron have been scaled 
down with respect to the relative solar abundance ratio to oxygen. The verti- 
cal line indicates the best fit temperature (2.2 X 10 6 K) in the isothermal case 
for the absorbing gas toward 4U 1820-303. 

density along the line of sight, are less affected by the differ- 
ent models and the charge exchanges between the solar wind 
and the local neutral gas, and, therefore, should provide a re- 
liable measurement of the iron abundance in the hot ISM. 

In this Letter, we present the first detection of the Fe xvn 
absorption line at ~ 15 A from the hot ISM toward 4U 1820- 
303. Comparing the column density of Fe xvn with that 
of Ne ix, and with those of O vn and O vm that we previ- 
ously measured, we have derived a relative abundance ratio 
of Fe/Ne and then inferred the iron depletion level in the hot 
gas. Throughout the Lett e r, we adopt the solar abundances 
from lAnders & Grevessel (1 19891) . 3 and quote the errors at 
90% confidence levels for single varying parameters until oth- 
erwise specified. 

2. SOURCE, OBSERVATIONS, AND THE EXISTING RESULTS 

4U 1820-303 is a bright low mass X-ray binary residing 
in the globular cluster NGC 6624 (Galactic coordinates l,b = 
2.°7 9,7.°91) and its distanc e has been determined as 7.6 ±0.4 
kpc dKuulkers et al. 12 003). meaning that it is very close to the 
Galactic center and is located ~ 1 kpc above the disk plane. 

The Chandra X-Ray Observatory has observed this source 
three times with different instrumental configurations. The 
observation log, data reduction and analysis procedures have 
been described in detail in Yao & Wang (2006, hereafter Pa- 
per I), and here, we summarize absorption line detections and 
relevant absorption line diagnostic results. 

We have detected highly ionized O vn, O vm, and Ne ix Ka, 
and O vn K/3 absorption lines, which are produced in the hot 
ISM rather than in the circumstellar gas associated with the 
binary system (Paper I; see also Futamoto et al. 2004; Juett 
et al. 2006). Modeling these lines with our absorption line 
model, absline (lYao & Wang 1120051) . we have constrained dis- 
persion velocity of the hot gas [v?, = 255(165,369) km s" 1 ], 
column densities of O vn, O vm, and Ne ix. We have also 
obtained the abundance ratio of Ne/O in the hot gas, which 
is consistent with the solar value. For a gas at tempera- 
ture about 10 6 — 10 7 K, the population of each abundant iron 
ion contained in the gas, e.g., Fe xvi, Fe xvn, and Fe xvm, 

3 The "real" solar abundances are still under debate (e.g., Asplund et 
al. 2005; Antia & Basu 2006), so we still use the old values. 
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A (A) 



Fe XVII 
Fe XVn 
Fe XVn 
Fe xvm 
Fe XVIII 
Fe XVIII 



2.v 2 2p 6 (' S)-2s 2 2p 5 4d( 1 P°) 

2i 2 2p 6 ('S)-2i 2 2p 5 3J( 1 P°) 

2s 2 2p 6 ( 1 S)-2s 2 2p 5 3d( i D°) 

2s 2 2p 5 ( 2 P°)-2s 2 2p 4 ( l S)3d( 2 D) 

2s 2 2p 5 ( 2 P°)-2s 2 2p 4 ( 1 D)3d( 2 P) 

2s 2 2p 5 ( 2 P°)-2s 2 2p 4 3d( 2 D) 



12.123 
15.015 
15.262 
14.121 
14.203 
14.361 



0.53 
2.31 
0.63 
0.90 
0.57 
0.93 



is distributed in a narrow temperature range (Fig. [TJ; there- 
fore, a well confined gas temperature or its distribution is 
crucial for inferring the total iron in the hot gas. The de- 
tection of multiple absorption lines in this sight line enables 
us to obtain such a constraint. For instance, if the interven- 
ing gas is in the collisional i onization equilibrium (CIE) state 
( Arnaud & Raymond! 1 1992h and isothermal, its temperature 
can be determined as T = 2.2 ± 0.3 x 10 6 K (Paper I). 

3. Fe XVII ABSORPTION LINE AND IRON ABUNDANCE IN THE 
HOT ISM 

We searched for the ionized iron absorption lines at the cor- 
responding rest frame wavelengths in the wavelength range 
between 2 and 25 A in the spectrum obtained in Paper I. The 
Fe xvi absorption lines are expected to be very weak (the os- 
cillation strength fii< 10~ 6 ), and are not considered further 
in this work. Table Q] lists the strong lines (/jj > 0.5) of ions 
Fe xvn and Fe xvm. We only detected a significant Fe xvn 
absorption line at 15.02A (Fig. |2), and did not see any clear 
sign for the other lines listed in Table [T] These detection re- 
sults are not surprising. The constrained hot gas temperature 
favors the Fe xvn population, and the transition of the 15.02 
A line is strongest (Fig. [TJ Table [TJ. Therefore, the Fe xvn 
absorption line at ~ 15.02 A is expected to be at least 4 times 
stronger (in terms of equivalent width [EW]) than the others. 

We use different models to characterize the Fe xvn absorp- 
tion line at 15.02A. The negative Gaussian model gives the 
line centroid as 15.008(14.999, 15.018) A or 132(-50, 314) 
km s , line width as a < 280 km s -1 or Vb < 396 km s , and 
its EW as 5.1(2.9, 7.3) mA. 

We then fit the line with our absline model, which adopts 
the Voigt function as line profile and allows a joint analysis 
of multiple absorption lines (see Yao & Wang 2005; Paper I; 
Wang et al. 2005 for further discussion). The fit is as good as 
with a Gaussian model, and the obtained line position is also 
identical. Connecting the v/, with those of O vn, O vm, and 
Ne rx lines (since the non-thermal broadening dominates, we 
therefore ignore the tiny differences of the thermal broadening 
in different elements), we obtain the column density of the 
Fe xvn as log|W F exvii(cirr 2 )] = 15.0(14.7, 15.2). 4 

Following the procedure we established in Paper I, we 
probe the abundance ratio of Fe/Ne in the hot gas. Since 
neon is a noble element, it is very unlikely depleted into dust 
grains. Therefore we take it as the reference element. In fact, 
we have obtained a Ne/O ratio that is consistent with the solar 
value (Paper I), and the following inferred Fe/Ne is essentially 
the same as Fe/O in units of solar value. Assuming that the 
absorbing gas is in a CIE state and isothermal, we jointly an- 

4 This column density, together with the dispersion velocity vj, constrained 
in jointly analyzing oxygen and neon lines (§2), can reproduce the EW mea- 
sured with Gaussian model. 
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alyze the Fe xvn line with the O vn, O vm, and Ne ix Ka lines 
and the O vn K/3 line, requiring the common absorbing gas to 
be of the same temperature. We fix the neon abundance at the 
solar value, and let the abundances of oxygen and iron be free 
parameters. In this way we constrain the hydrogen column 
density to A? H = 7.9(5.0, 10.2) x 10 19 cm" 2 , the temperature to 
T = 2.2(1.9,2.5) x 10 6 K, and the abundance ratio to O/Ne = 
0.9(0.5, 1.3) solar for the absorbing hot gas, which are iden- 
tical to those reported in Paper I. In addition, we obtain the 
abundance ratio of Fe/Ne as 0.8(0.4, 2.1) solar. Considering 
the dependence of the Fe xvn population on T (Fig. [T), we 
calculate the confidence contours of Fe/Ne versus T, which is 
presented in Figure^. 

Next, we investigate the effects on the inferred Fe/Ne ratio 
if the above isothermal assumption of the intervening gas is 
relaxed. Since the absorption samples almost the entire Galac- 
tic disk from the Sun into the Galactic bulge, it is possible 
that the hot gas consists of different temperature components. 
Here we examine two simple temperature distributions. In 
each case, we first interpolate the ionization fractions at dif- 
ferent temperatures, assuming the gas to be in CIE state, and 
then calculate the column density for each ion. To get a bet- 
ter constraint, we also add the undetected Ka Ne x (12.134 
A) absorption line in the joint fit. This line, except for the 
detected line absorptions from oxygen, neon, and iron ions, 
is the next most expected one from a different ion to be ob- 
served in the spectrum with high counting statistic because of 
its anticipated large column density (Fig. [TJ and large transi- 
tion coefficient (/ij = 0.416), and is particular useful for con- 
straining the upper boundary to the gas temperature. In the 
first case, we assume that the hot gas temperature distribution 
follows a logarithmic Gaussian form, as a natural extension of 
the isothermal single temperature case, 



dNn(T) oc exp 



-( logr- log7b) 2 



dlogCT), 



(1) 



2(Oi og 7-) 2 

where the mean temperature Tq is equivalent to T in the 
isothermal case, and o"i og r is the dispersion of log To. Un- 
der this assumption, we obtain 7q = 2.0(1.8,2.4) x 10 6 K and 
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FIG. 2. — Fe XVII absorption line at 15.02A in the Chandra MEG-LEG 
combined spectrum of 4U 1820-303, modeled with the absline model (his- 
togram). The dip at 15.26 A indicates the expected amount of absorptions for 
the second most significant Fe XVII line (Table[TJ. The dotted line indicates 
the smooth continuum level around the line. The model has been convolved 
with the instrumental response, and the binsize is 12.5 mA. The apparent 
bump on the blue-side of the line may be due to the imperfect calibration 
near the hot pixel. For the detail information of the co-added spectrum and 
the continuum fit, please refer to Paper I. 
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FIG. 3.— The 68%, 90%, and 99% confidence contours of (a) the abun- 
dance ratio of Fe/Ne (in units of solar value) vs. gas temperature T, (b) Fe/Ne 
vs. temperature dispersion, and (c) Fe/Ne vs. power-law index 7 for the 
isothermal case, log-Gaussian temperature distribution, and the power-law 
temperature distribution, respectively. 

cio g r < 0.15, and the abundance ratio Fe/Ne = 0.8(0.4, 2.3) 
solar. In Paper I, we have demonstrated that since an isother- 
mal absorbing plasma is adequate to describe the observation, 
the additional free parameter cr log T cannot be fully constrained 
in the spectral fitting, and there is an apparent correlation be- 
tween To and <r\ og T (see Fig. 4 in Paper I) due to the large 
value of A^ovm- But the inferred Fe/Ne is insensitive to differ- 
ent (Ti og r values, as illustrated in the confidence contours of 
the Fe/Ne ratio versus (Ti og r (see Fig.[3j3). 

In the second case, we assume that the hot gas temperature 
distribution follows a power-law (PL) form, 



dN H (T) : 



AT H q(7+ 1) 
To 



(T/ToVdT. 



(2) 



This simple characterization of the temperature distribution 
can be deri ved, for example, n aturally from an exponential 
disk model dYao & Wang 1120071) . where Nno is the total hy- 
drogen column density along the sight line, To is the Galac- 
tic mid-plane temperature, and the PL index 7 is the ra- 
tio of the temperature to density scale heights. Our joint 
analysis gives N H = 7.5(5.3,10.0) x 10 19 cm" 2 , 7 > 2, and 
T = 2.4(2.1,3.4) x 10 6 K. Again, although the extra free pa- 
rameter 7 can vary in a large range, the constrained abundance 
ratio Fe/Ne=0.9(0.4, 2.0) solar is insensitive to 7. Fig. [3J: 
shows the confidence contours of the Fe/Ne ratio versus the 
PL index 7. 

4. DISCUSSION 

We detect a significant Fe xvn absorption line at ~ 15.02 
A in the Chandra spectrum of 4U 1820-303. A joint anal- 
ysis of this line with the detected highly ionized oxygen 
and neon lines, all interstellar in origin and observed in the 
same spectrum, gives the abundance ratio of Fe/Ne in the hot 
ISM component, which, although with large errors, is con- 
sistent with the solar value. In addition, this result appears 
to be unaffected by the different gas temperature distributions 
adopted. We conclude that there is no evidence for substan- 
tial depletion of iron into dust grains in the hot ISM. This 
is in contrast to cooler phases of t he ISM where iron is usu- 
ally found to be heavily depleted (Sembach & Savage Hl996l : 
lJuett et aL~ll2006l) . Grain cores containing iron oxides are gen- 
erally rather resili ent and it is quite difficult to liberate the iron 
from these cores (Sembach & Savage 1996; Frisch & Slavin 
120031) . This solar value of Fe/Ne ratio, if confirmed, thus sug- 
gests that likely all of the dust in this very hot ISM phase (T 
> 10 6 K) has been destroyed by frequent and/or severe shocks 
during the dust grain processing in the ISM. 

Dust grains pre-existing in the ISM or formed in supernova 
ejecta can be destroyed by their generated forward and re- 
verse shocks and subsequently, in heated hot gas. By study- 
ing the iron abundance via the far-UV Fe 11 absorption lines, 
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ISembach & Savage I (1 19961) find that while more than 99% of 
iron is depleted into dust grains in the cold disk of the Galaxy, 
the iron depletion is ~ 80% in warm clouds of the Galac- 
tic halo. They attribute this difference to the dust grain dis- 
ruption by the supernova (SN) shocks when circulating the 
grains between the Galactic disk and halo. If the hot ISM is 
believed to be heated from the cool ISM, it should have expe- 
rienced shocks much more frequently and/or much more vio- 
lently than the warm halo clouds. Therefore it is natural to ex- 
pect that many more grains, even including the resilient iron- 
rich cores, could have been destroyed in such harsh environ- 
ments. More recently, Strickl and et al. I (120041) have obtained 
an abundance ratio of Fe/O from the diffuse extraplanar halo 
emission of many disk galaxies like our own, which is ~ 40% 
solar. This result, although subject to different thermal plasma 
models adopted and the different thermal properties assumed 
for the emitting hot gas in the data analysis, clearly rules out 
the depletion pattern found by Sembach & Savage] (119961) in 
the cold and warm gas of the Galactic disk and halo, further 
supporting the scenario that more iron has been released back 
to gas phase. 

The above interpretation may not be entirely unique and 
in some cases is subject to systematic effects introduced by 
line-of-sight variations and intrinsic X-ray source properties. 
Claims like an overabundance of heavier elements relative to 
oxygen in neutral matter towards the Galactic center direc- 
tion presumably caused by a significant contr ibution of Type 
la SNe in the Galactic bulge (Ueda et al. 2005) would elevate 
the contribution o f iron in the hot ga s phase as well. In this re- 
spect the result by lJuett et al. I (120061) that iron is significantly 
depleted in the cool phase in line of sight towards 4U 1820- 
303 and that Ne/O is close to the solar value is an important in- 
dicator that conditions are not so unusual. Abnormal heavy el- 
ement abundances have been observed in microquasars (e.g., 
Cyg X- 1 : Schulz et al. 200 1 , Juett et al. 2004; GRS 1 9 1 5+ 1 05 : 
Lee et al. 2002, GX 339-4: Miller et al. 2004). In all these 
systems, the observed overabundance pattern of the heavy 
elements relative to oxygen can also be interpreted as the 
different photoionization struc tures of these elements in the 
intrinsic absorbing material ( Sch ulz et al. 1 120021 : iLee et al.~l 
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120021) . iFutamoto et ail (2004) and Juett et al. (2004, 2006) 
already argued successfully that this is likely not the case for 
4U 1820-303. 

As a final note, we point out that the absorption path length 
toward 4U 1820-303 passes through the Galactic bulge re- 
gion where the soft X-ray bac kground emission in th e 0.75 
keV band is greatly enhanced (Snowd en et al. l!T997t) . This 
enhancement indicates that the emitting plasma is either of 
a high temperature or of a dense emitting region, or both. 
Recently, we have obtained an exponential scale height of 
~ 2 kpc and the mid-plane density of ~ 2.4 x 10~ 3 cm -3 
for the hot gas toward Mrk 421 (Galactic coordinates l,b = 
179.°83,65.°03) in a joint-a nalysis of absorptio n and emis- 
sion data at the same time dYao & Wang Il200~7l) . This char- 
acterization, if typical for the general hot ISM, only accounts 
for ~ 30% of the observed absorption toward 4U 1820-303, 
meaning that a large portion of the absorption originates from 
the Galactic bulge. Therefore, the gas phase abundance ratio 
of Fe/Ne presented in this letter could be biased because of the 
remarkable absorption contribution from the bulge, where the 
metal abundance pattern may not be as the same as that in the 
overall ISM. Nevertheless, we present here a feasible way to 
infer the gas phase iron abundance in the hot ISM that poten- 
tially affects our understanding of the cooling/heating process 
in the ISM in general. To obtain a global picture of the gas 
phase iron abundance in the hot gas, high quality absorption 
data along other sight lines that are away from the Galactic 
bulge region are therefore required. 
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